Broken time-reversal symmetry in strongly correlated ladder structures 
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We provide, for the first time, in a doped strongly correlated system (two- leg ladder), a controlled 
theoretical demonstration of the existence of a state in which long-range ordered orbital currents 
are arranged in a staggered pattern, coexisting with a charge density wave. The method used 
is the highly accurate density matrix renormalization group technique. This brings us closer to 
recent proposals that this order is realized in the enigmatic pseudogap phase of the cuprate high 
temperature superconductors. 



The circulating current phases in correlated electron 
systems, also called orbital antiferromagnets (OAF), 
were first considered in the context of excitonic insulators 
EJ, but then discarded in favor of more conventional or- 
der. After the discovery of the cuprate high temperature 
superconductors they were rediscovered H Q and called 
staggered flux (SF) phases Q. Many of their properties 
were discussed, but were forgotten again in the absence 
of experimental vindication. The discoveryof an unusual 
and robust regime called the pseudogap [4| in these su- 
perconductors has changed the picture once more. The 
pseudogap mimics the momentum dependence of the su- 
perconducting gap, but the state itself is not supercon- 
ducting. In this context, two recent developments have 
taken place: (1) attempts have been made to explain this 
regime in terms of fluctuations of SF order and (2) a 
proposal has been made that it is not fluctuations, but 
a true broken symmetry that is the origin of the pseudo- 
gap H ■ This ordered state was called the singlet d- 
density wave (DDW) following Ref. @ and is the same 
as the OAF and SF phases. In this Letter we adopt 
the density wave (DW) terminology, as it can describe 
large classes of order parameters with orbital angular mo- 
mentum. The label d stands for angular momentum 2, 
as in atomic physics. The conventional charge density 
wave (CDW) in which charge is modulated in space is its 
s-wave counterpart with angular momentum zero. The 
triplet s-wave density wave is what is commonly called 
a spin density wave (SDW). Another breakdown of time 
reversal symmetry in which the circulating currents obey 
translational symmetry, as opposed to DDW, has also 
been pointed out and has been argued to be responsible 
for the pseudogap phase H 0] . 

Although much indirect experimental evidence of 
DDW can be argued to exist, a direct observation of 
DDW would be Bragg reflection of neutrons carrying 
magnetic moments from the staggered arrangement, on 
the scale of a few A, of circulating currents. Recent 
neutron scatte ring experiments have, however, been con- 
troversial E2> El El E~l and more precise and well- 
characterized experiments are underway to establish this 



FIG. 1: Circulating plaquette currents on a ladder, charac- 
teristic of the DDW phase. 



order. Thus, theoretical exploration of microscopic mod- 
els of correlated electronic systems with controlled meth- 
ods has acquired urgency. We, therefore, study the sim- 
plest geometrical structure in the form of a two-leg ladder 
[15| that can support staggered orbital currents, as shown 
in Fig. Q] 

Previous studies of DDW order in two-leg ladders have 
used weak-cou pling bosonization /renormalization group 
(RG) analyses El El El El Efl El 1110, density- 
matrix renormalization group (DMRG) [2J| analysis of 
the i-J-model [23| and a half-filled Hubbard-like model 
|25). or exact diagonalization [2(J of the t-if- J-model. At 
half-filling, models with long-range ordered currents have 
been found both for spinless E3 and spinful [H 0, |H 
lifij l fermions. In contrast, for doped ladders, power-law 
decay has been found for spinless 0, and spinful 
El [3 IH cases. For the t-J ladder, only short-range 
order was found [l^, and the study of the t-t'-J model 
[2f| yielded similar results. 

The approach used in the present work is the accurate 
DMRG method that can be used for arbitrary interaction 
strength, unlike the weak-coupling bosonization/RG ap- 
proaches. The results of our calculations are striking. 
Although common t- J-type models do not exhibit long- 
ranged DDW order, a separate class of repulsive Hamil- 
tonians show robust long range DDW order even in the 
presence of substantial doping. These have their histor- 
ical origin in a half-filled SO(5) invariant model on a 
ladder [27j. At precise half- filling, it was shown to ex- 
hibit DDW in its phase diagram (2^, |2|j . We shall show 
that SO(5) invariance is irrelevant by considering cou- 
pling constants very far from the "SO (5) parameters" 
and by substantially doping this model. The real rea- 
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son for success with this class of models is that it strad- 
dles CDW and DSC- like (d 2 .2_j,2-superconductor) states 
(more precisely rung-singlet states), resulting in a local 
kinetic exchange between them. This is much like the ac- 
tual situation in the cuprates in which the DDW phase is 
an intermediate regime between a multi plic i ty o f complex 
charge ordered states and DSC d |H M 13 . 

We label the site of a ladder by i = (r, I), where r = 
1, . . . , L is the rung index and I — 1,2 is the leg index. 
The current operator between any two sites i and j, J7y, 
is 

(7 

where c\ is the creation operator of a fermion with spin 
a at site i. We set the lattice spacing to unity. 

There are at least two convenient ways of probing 
DDW order. One approach is to measure the equal- 
time rung-rung current correlation function in the ground 
state 

C(r, R) = (j Iuns (R + r/2)j xung (R - r/2)}, (2) 

where jmngM = c7(r,i),(r,2V I n order to minimize the ef- 
fect of the boundaries of a finite ladder, we should choose 
R to be the location of the central rung, and we shall de- 
note this correlation function as C(r) 

An alternative approach is to break the time reversal 
symmetry explicitly by applying a source —hj lung (l) on 
one end of the ladder and measure the current induced 
in the sample. The source term for DDW is necessar- 
ily complex and a complex DMRG program is needed, 
which is more demanding on memory and computer time. 
Nonetheless, we have used both methods for cross checks 
for every single case discussed in the present Letter. Of 
course, the sensitivity with respect to the magnitude of 
the source must be tested, and it must be made sure that 
it is not so large that we are probing some excited states. 
We check this issue very carefully in all our numerical 
work. For both methods, we use a finite size algorithm, 
which is more reliable, performing sweeps to reach con- 
vergence 0. 

To orient ourselves we shall begin with the two-leg t-J 
ladder, which is defined by the Hamiltonian: 

Hu = -t K^P+Lcl + J^fSi.Sj-^), 

(U>,<r <ij> 

(3) 

where n\ is the total fermion occupation number at the 
site i, and Si is the spin operator at the corresponding 
site. V projects out doubly occupied sites and de- 
notes pairs of nearest neighbor sites. The t-J model is 
the simplest model that captures some aspects of the 
high temperature superconductors. Removal of electrons 
from the system is quantified by the doping parameter 
(5=1 — (rii). In actual experiments, DSC is observed 
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FIG. 2: The correlation function C(r) of the t-J-V-V' model. 
The parameters are J/t = 0.4, V/t = 3, V /t = 1, and 8 = 
0.1. The size of the ladder is 120x2. We kept 800 states and 
performed 42 sweeps. The vertical scale is chosen in units 
where t =1. 

in the range S = 0.05 — 0.25; at small values of 5 lies 
the pseudogap regime. In this model, the DDW cor- 
relations decay exponentially with a correlation length 
£ = 3 — 4 . Additional next-nearest neighbor hopping, 
augmenting the model, may be supposed to suppress the 
competing CDW and DSC phases and thus reveal DDW 
order, but our calculations and those of Ref. |2(| do not 
support this idea. 

A more interesting model is the t-J-V-V model, which 
is a t-J model augmented by Coulomb repulsion terms 
V and V between nearest and next nearest neighbors: 

HtJW = Hu + V^2 n i n i + V ' ni "j' ( 4 ) 

(U) «y» 

where ((i, j)) denotes next nearest neighbor pairs of sites. 
In order to reduce the effect of open boundaries, a chem- 
ical potential term (V + V')m is added to the bound- 
ary sites. A typical parameter set, given by J/t = 0.4, 
V/t = 3, V'/t = 1, and S = 0.1, yields the current 
correlation function C(r) shown in Fig. We observe 
what appears to be a bubble of DDW extending up to 
20 rungs, and then a sinusoidally modulated exponen- 
tial decay with £ « 3. Thus, although we observe DDW 
over a moderately long range, there is no macroscopic 
order. We have probed this model by the second method 
in which we induce a current by a source at the edge. As 
shown in Fig. |21 the long range correlations that are ob- 
served at the first infinite-size step disappear after a few 
sweeps and converge after six sweeps to an exponential 
decay with a correlation length ( « 10. The moderately 
large bubble- like nature of C(r) and the different values 
of £ in these two methods are strong reasons for suspect- 
ing proximity to a first order transition to the ordered 
DDW phase, as the boundary effects seem to nucleate 
this phase. 
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FIG. 3: Absolute value of the current summed around a pla- 
quette, |0'piaq(j"))|, as induced by an edge current, h — O.Oli, 
as a function of the distance r from the edge of a 60 x 2 lad- 
der in the t-J-V-V model, with the parameters described in 
the text. Note the convergence as a function of sweeps. The 
number of states retained were 400 and the results were found 
to scale with current at this level. 



Finally, we shall consider a different class of Hamilto- 
nians, in which a pair of electrons across a rung is given 
an internal structure, much like a molecule. This is an 
interesting way of generating a set of low energy Hamil- 
tonians |27j that are defined by 



[J±S rA ■ S r , 2 + V±{n r , 1 - 1)K, 2 - 1)] . (5) 



Note that there are no longer any projection operators 
V, as in the previous examples and the problem can be 
treated for arbitrary interaction strength. In the half- 
filled case this Hamiltonian has a precise SO (5) symme- 
try Eifl w hen J ± = 4(17 + VI). It also exhibits DDW 
order |2ll Ej| . But as soon as the system is doped, or 
the parameters are no longer finely tuned, there is no 
SO (5) symmetry. The weak coupling phase diagram at 
half- filling obtained from bosonization/RG, as shown in 
Fig- EI gives us some guidance as to where to look in our 
DMRG calculations. This phase diagram is essentially 
identical to that of Ref. |22j, except that we also show 
the regime for U < 0. Other than the DDW and CDW, 
there are two relevant states that can be adiabatically 
continued to resonating valence bond states j^j °f the 
short-range variety 29] . These states are represented by 

M 



I rung-singlet) cx 









T/r. 



4 
3.5 

3 
2.5 

2 - 
1.5 - 

1 - 

0.5 - 

— 
-1 



SO(5): J,_ = 4(U + VJ 



U + 3/4 J ± = V ± 



Rung singlet 



Site 
singlet 





u/v ± 



FIG. 4: The weak coupling phase diagram at half filling for 
Htj ± uv ± from a bosonization calculation. The open circles 
correspond to the parameters in Fig.|S]and the open triangles 
to those in Fig. |S| 
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The DDW lies between the CDW and the rung-singlet 
phases. 

We have studied the Hamiltonian in Eq. JSJ for a 
range of parameters and find long-range DDW order in 
the doped model, which has nothing to do with SO (5) 
symmetry. Nonetheless, the DDW phase is situated be- 
tween the CDW and the rung-singlet phases, similar to 
the weak coupling bosonization results for half-filled lad- 
ders. As a typical example, we have shown in Fig. [S] our 
results for the rung current as induced by an edge cur- 
rent of tiny magnitude O.OOOlt. The parameters chosen 
were U = 0.25, t = V± = 1, J± = 0.8,1.1,1.5,1.7, and 
S = 0.04. As a response, we see robust long-ranged DDW 
order in the middle of this range of J± with stripe like 
features where pairs of holes reside; see, in particular, 
Fig. [5f c) , where we also plot the hole density, and the 
coexistence with stripe order is especially evident from 
the antiphase domain wall structure |3lj . The induced 
currents clearly alternate in sign and can be of order 
unity, in units of t, even though the source current is 
infinitesimally small. For ladders of lengths 100, 150, 
and 200, and for parameters of Fig. 02c), the current 
amplitudes are respectively 0.56, 0.53, and 0.53 respec- 
tively, consistent with long range order, though in a nu- 
merical calculation it is never possible to rule out a very 
slow decay. We have studied the d-wave pairing corre- 
lations, and find only extremely rapid decay in ladders 
that exhibit DDW long-range order. This is as expected 
from previous numerical work on the half-filled system 
[25J where it was shown that the superconducting corre- 
lations decay exponentially in the DDW phase. It is also 
in accord with the phase diagram in Fig. 0] which shows 
that the region with strong d-wave superconducting cor- 
relations, the "rung-singlet" phase, is distinct from the 
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FIG. 5: Rung current JrungfV ) as a function of the location of 
the rung r in a t-J±-U-V± model at 4% doping on a 100x2 
ladder, with parameters U = 0.25, t = V± = 1 and an edge 
current of 0.0001. The sequence of figures correspond to (a) 
Jx = 0.8, (b) J x = 1.1, (c) Jx = 1.5, and (d) J ± = 1.7. 
In (c), we show the profile of the hole density depicted as 
solid dots corresponding to the scale on the right. We kept 
up to 400 states and performed up to 8 sweeps. Note the vast 
differences in the scales of the current strengths. 

DDW region. More generally, bosonization leads to the 
prediction that the phase with strong DDW correlations 
has only short-ranged pairing correlations. For the case 
of Fig. Efc), we have also studied the low-lying excita- 
tions about the ground state. We find a robust spin gap, 
again in accord with general expectations. In Fig. El we 
show the results for 5 = 0.08 and for U = 0.5. For suffi- 
ciently strong doping, roughly between 10 to 20%, DDW 
is suppressed for these sets of couplings. We have con- 
firmed the existence of DDW on longer ladders, up to 200 
rungs, with current amplitudes decaying only slightly at 
the boundaries. It should be noted, however, that in 
two-leg ladders, there cannot be nodes in the excitation 
spectrum, as in two dimensions, which can, in principle, 
change the detailed nature of the doped phase. 

In summary, we have shown that there are repulsive 
microscopic models that exhibit DDW order at finite 
doping, providing added support for the identification of 
the pseudogap in the cuprates with this state. Our work 
also raises the real possibility that a n ap preciation of the 
complexity of many novel materials [32| may be impossi- 
ble without these remarkable broken symmetries, and it 
is important to search for such complex quantum order 
in an even wider class of Hamiltonians. 
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Gerhard-Hess prize of the Deutsche Forschungsgemein- 




5 10 15 20 25 30 35 40 45 50 5 10 15 20 25 30 35 40 45 50 



FIG. 6: Rung current j lmg ( r ) as a function of the location 
of the rung r in a t-J±-U-V± model on a 50x2 ladder, with 
parameters t = V± = 1, Jx = 1.2 and an edge current of 
0.0001. The sequence of figures correspond to (a) U = 0.25 
at 8% doping, and (b) U = 0.5 at 4% doping. We kept up to 
400 states and performed up to 8 sweeps. 
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